The entanglement dynamics of interacting qubits embedded in a spin environment 

with Dzyaloshinsky-Moriya term 
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We investigate the entanglement dynamics of two interacting qubits in a spin environment, which 
is described by an XY model with Dzyaloshinsky-Moriya (DM) interaction. The competing ef- 
fects of environmental noise and interqubit coupling on entanglement generation for various system 
parameters are studied. We find that the entanglement generation is suppressed remarkably in 
weak-coupling region at quantum critical point (QCP). However, the suppression of the entan- 
glement generation at QCP can be compensated both by increasing the DM interaction and by 
decreasing the anisotropy of the spin chain. Beyond the weak-coupling region, there exist resonance 
peaks of concurrence when the system-bath coupling equals to external magnetic field. We attribute 
the presence of resonance peaks to the flat band of the self-Hamiltonian. These peaks are highly 
sensitive to anisotropy parameter and DM interaction. 

PACS numbers: 03.65.Yz, 05.40.-a, 03.67.Mn, 75.10.Pq 



INTRODUCTION 



Since the discovery of quantum mechanics, quantum 
entanglement has played an important role in quantum 
information processing (QIP), such as superdense cod- 
ing teleportation [2] and quantum algorithms [3|. In 
recent years, there has been a growing interest in entan- 
glement dynamics by using the coherent manipulation in 
solid state systems. One of the most natural candidates 
is spin chain, which can be simulated by means of cold 
atoms in optical lattices [4] or by coupled microcavities 
5], and has been extensively studied in numerous works 
6-10]. Especially much attention has been paid to the 
relation between entanglement and quantum phase tran- 
sitions (QPTs) which are driven purel y by quantum fluc- 
tuation occurring at zero temperature jlll-[l3|. The QPT 
is related to a dramatic change in the ground-state prop- 
erties of the system as an external non-temperature pa- 
rameter varies across the transition point. Consequently, 
the dramatic change in the structure of the ground state 
should result in a great difference between the quantum 
correlation on both sides of the quantum critical point 
(QCP). Naturally, the entanglement should be able to 
characterize the change. For example, Osterloh et al. \J} 
has proven that the derivative of the nearest-neighbor 
entanglement diverges at QCP. Due to the dynamical ul- 
trasensitivity of the induced quantum critical system [ 14] , 
quantum entanglement can be treated as a tool [iBI^Ilq 
to characterize QPTs. 

It's worth noting that the systems considered are 
closed, i.e., isolated systems have no interaction with 
their external environment. However, real physical sys- 
tems are never isolated, since the coupling between sys- 
tem and the surrounding environment is inevitable. The 
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quantum dynamics of physical systems is always compli- 
cated by their coupling to many 'environmental' modes. 
The dominant environmental effects are localized modes 
at low temperature, which are usually described by spin- 
bath model. An interesting phenomenon is that the cou- 
pling process between system and bath shows a duality of 
influence on quantum system. On one hand, the coupling 
can assist people to achieve some tasks in QIP [l^, • 
One of the focuses is the induced entanglement between 
the two noninter acting qubits. For instance, Yi et al. 
[19] showed that the entanglement changed dramatically 
along the line of critical points of spin bath. Subse- 
quently, an enhanced effect of induced entanglement near 
the critical point was demonstrated ^]. These investi- 
gations exhibit a new perspective to engineer protocols 
for entanglement generation. On the other hand, the 
coupling between system and bath can play a role as 
decoherence. It can transfer a pure ensemble of qubits 
to a mixture of classical ones [21] and lead to asymp- 
totical disappearance of system entanglement. In some 
cases, the entanglement will vanish even in finite time 
[22] . Great efforts have been devoted to the study of the 
decoherence caused by the spin bath [23, 24]. One of the 
most common models is anisotropic XY model, which en- 
compasses two well-known spin models, i.e., Ising chain 
and the XX (isotropic XY) chain in a transverse field. 
The XY model holds an advantage that it can be exactly 
solved by mapping to a spinless fermionic model. Such 
solvability provides a playground for testing the physi- 
cal ideas [25]. For one-qubit case, Quan et al. [14] first 
proved that the decay of the Loschmidt echo (LE) was 
enhanced by the QPT of the Ising bath. With that, this 
consideration was extended to the two-qubit case. In the 
transverse Ising model. Sun et al. [26] showed that the 
concurrence decayed exponentially with fourth power of 
time in the vicinity of the critical point of spin bath. 
For an XY spin chain with Dzyaloshinsky-Moriya (DM) 
interaction, Cheng et al. [27] found that decay of deco- 
herence factor was sensitive to the DM interaction, espe- 
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cially in the strong-coupling region. Also, the three-qubit 
case was studied by some works [28|, [29| . The previous 
works have extensively studied the decoherence process 
of initially entangled state of noninteracting qubits |26|- 

In this paper, we investigate the influence of the XY 
spin bath with DM interaction term for two interacting 
qubits. The paper is organized as follows. In Sec. [Ill we 
introduce the model Hamiltonian describing two-qubit 
system coupled to an XY spin chain with DM interaction, 
and derive an analytic formula for the exact solution of 
entanglement. In Sec. IIIIl regarding the spin bath as a 
decoherent environment, we exhibit the competition be- 
tween system-bath and interqubit interaction in the en- 
tanglement dynamics. We analyze the competing effects 
of environmental noise and two-qubit interaction for var- 
ious system parameters, especially at the QCP of spin 
bath. We observe resonance peaks emerging for some 
specific parameters. Finally, we give a summary of our 
results in Sec. HVl 

Environment 
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FIG. 1: (Color online) The schematic diagram of two exter- 
nal interacting spin qubits symmetrically coupled to the envi- 
ronment described as an XY spin chain with Dzyaloshinsky- 
Moriya term. 
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where ctj {a = x, z) is the Pauli matrix at site j. 
The parameters 7 and A characterize the anisotropy of 
the self-Hamiltonian and the intensity of the magnetic 
field applied along the z axis, respectively. Here, such 
self-Hamiltonian comprises DM term, which is an 
antisymmetric spin coupling. The DM interaction often 
arises from a mixture of superexchange and spin-orbit 
coupling in low-dimensional magnetic materials |32| . 
For analytical solvability, we assume that the D vector 
is imposed along the z direction., i.e., D = Dz. The 
strength of interaction between the system qubits is given 
by J(> 0), and the coupling between the system qubits 
and the surrounding spin chain is denoted by g. 

Since [Hab^ Hj] = 0, Hamiltonian ([1]) can be rewritten 

as 



H = Y,eM{4>^\®H^E'' 



(5) 



where (/i = 1, ... ,4) is the /ith eigenstate of Hab 
with eigenvalue e^. Under these bases, the coupling Hj 
between the qubits and spin chain exerts an extra mag- 
netic field on the spin bath H^, and then gives rise 

to an effective qubit-dressed Hamiltonian H'^^ with mag- 
It is easy to find that 



netic field = A 



6 



0, ^4 



(6) 



II. MODEL AND SOLUTION 



The Hamiltonian H^'^ can be diagonalized by following 
the standard procedures. The Jordan- Wigner transfor- 
mation [H, [13 



The quantum system we consider consists of two inter- 
acting qubits and coupled to a spin environment 
E. The environment is composed by N spin- 1/2 par- 
ticles with periodic boundary condition, which can be 
described by one-dimensional XY spin chain with DM 
interaction [30]. The schematic diagram of the quan- 
tum open system is shown in Fig. [H The corresponding 
Hamiltonian is given by 



H = + Hab 



Ht 



(1) 



cr+ exp 



CF. exp 
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ct| = 1 - 2f]f^, 



(7) 



maps spin chain to one-dimensional spinless fermionic 
model with creation and annihilation operators as fol- 
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il-2iD)fjf,+^ + {l + 2iD)fj^Jj 

+ 7 {fjfli + h+^f^) + (l - 2/j/,) 
(1 - 2iZ))/t,/i + (1 + 2iD)flfN 
{fyi + exp(i^) + (l - 2fy, 



(8) 



An extra phase (j) = 7r(^^-,^ fjfj + 1) appears on the 
chain boundary due to phase accumulation of the Jordan- 
Wigner transformation. 

Next discrete Fourier transformation is introduced to 
convert the fermionic operators from real space to mo- 
mentum space by defining 



dk 



(10) 



with the discrete momentums as 
(2n7r + ^) _ TV TV 

After that, the Hamiltonian becomes 

H^^ = ^(2cos/c-2A^ -4Dsinfe)44 

k 

-^i-fsink (dlj^dl ^ d-kdk^ -^NX^. (11) 

The diagonalized form is achieved by Bogoliubov trans- 
formation which defines quasiparticle creation (annihila- 
tion) operator {^k,x^) as 
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with the angle O^^^ defined by 

.A, -J sink 



tan 0. 



Xj^ — cos k 



(12) 



(13) 



Consequently, the dressed self-Hamiltonian is unitarily 
equivalent to such diagonal form 



(14) 



where 



{k)=2^J (A^ - cos kf ^{j sin kf - 4D sin k. (15) 



The ground state |G)a^ has no quasiparticle for arbi- 
trary /c, i.e., bk^x^\G)x^ = 0. Due to the relation 

5fe,A, = cos(^^/2) dk-isin{0l^/2)dl^, |G)a, can be 
written as 

|G)a, = n (cos^|0)fc|0)_fe +isin^|l)fe|l)_fc J , 

k>0 \ ) 

(16) 

where |0)/c and |l)/c are the vacuum and single excitation 
of the /cth mode d/e, respectively. With these analyti- 
cal expressions, we can straightforwardly derive the time 
evolution of arbitrary initial state and obtain the reduced 
density matrix of the two-qubit system, and then exam- 
ine the effect of the environment. Suppose that at time 
t = the qubits are completely disentangled from the en- 
vironment, i.e., the global system wave function is given 

by 



i*(o)) = \m)AB ® \m)E. 



(17) 



Clearly, governed by the Hamiltonian ([T]), the state 
at time t is given by \^(t)) = /7(t)|^(0)), where U{t) = 
ex.p{—iHt) is the evolution operator of the composite sys- 
tem. The reduced density matrix of the two-qubit system 
is obtained 



-iHt 
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where 



^>^.mO))AB, 



(19) 



and the decoherence factors are 



i^M-(t) =£;((/>(0)|exp(ii/^-t)exp(-ii/^^t)|(/)(0))£;. (20) 

We assume that the two qubits in AB system ini- 
tially stem from a separable state, i.e., |0(O))ab = 
(|0) + |1))/V2 (g) (|0) + |1))/V2, where |0) and |1) denote 
the spin up and down, respectively. The initial state of 
the environment is supposed as the ground state of i?^, 
i.e., \(I){0))e = \G)x. By tedious calculation [13, |35|, we 
have 

FM.(t) = A(G|e^^-^^e-^^-'^^|G)A 



= n^osO^^ cose^e*(^^'-^'«^)^cos(6^^ -6^) 

k>0 

\Au /^Aiy\ 



sin sin e-'^"^" ""^ cos(efc^ - 6^ 
+ sinefe'' cose^''e-*("^+"^)*sin(efe'' -e^'') 



cos sin e^e^("'=''+"^)* sin(e^ - e^ ), 



(21) 
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where 0, 



-0^)/2 is the angle difference between 



the normal mode dressed by the system-environment in- 
teraction and the purely environment. 

From Eq. (p!8|) . in the bases spanned by 
{|00), |01), |10), |11)}, the reduced density matrix of the 
two-qubit system is obtained 



pAB{t) = - 



'^Fi2{t) 

1 

^2*3 (t) 

■^F,\{t) 



'^Fi3{t) 
F23{t) 
1 

-^F,\{t) 



Fl4(t) 

e—F2A(t) 

iJt _ , . 

e— F34(t) 



1 

(22) 

Some matrix elements can be further simplified for the 
choice of A^. For the case of Eq. (|6]), we have 0^^ = 
6^^ = 0, and also it is obvious that such relations hold: 
Fi2 = Fi3, F23 = 1, F24 = F34. So there are three 
independent decoherence factors. 

To investigate the entanglement dynamics of two 
qubits surrounding spin bath, rather than decoherence 
factor employed in Refs. [23, [35|, we utilize the con- 
currence directly [36^, an entanglement measure for any 
bipartite system that relates to the two-site reduced den- 
sity matrix p. The concurrence for two qubits is defined 
as 



(23) 



C{t) = max{0, ^/uJl — yJZJ^ - 



where uji (i = 1, . . . , 4) are the eigenvalues in decreasing 
order of the auxiliary matrix C = pioy (8) (7y)p*{(7y cr^). 
Here p* denotes the complex conjugation of p in the stan- 
dard bases and Gy is the Pauli matrix. The concurrence 
varies from C = for a separable state (zero entangle- 
ment) to C = 1 for a maximally entangled state. 



III. ENTANGLEMENT DYNAMICS OF THE 
TWO-QUBIT SYSTEM 

Before we consider the competition of system-bath cou- 
pling and two-qubit interaction, we first emphasize the 
following two limiting cases. On one hand, when there is 
no coupling between system and environment, i.e., ^ = 0, 
we can get uji = sin^ Jt/2, lj2 = ^3 = ^4 = 0, and then 
C{t) = I sin Jt/2|. It means the interaction between the 
qubits can generate an oscillating entanglement with pe- 
riod of 27r/ J. At t = (2m + l)7r/ J (m = 0, 1, 2, ...), the 
state of the system reaches the maximum entanglement. 
On the other hand, when there is no interaction between 
the two qubits, the entanglement dynamics of the qubit 
system returns to the process of entanglement induced by 
the spin bath pJi]. The entanglement changes dramat- 
ically along the line of critical points of the spin bath. 
Now we investigate the infiuence of the spin bath dur- 
ing the process of entanglement generation. In Fig. [2l 
concurrence is plotted as a function of magnetic field A 
and time t with J = 2,7 = l,D = 0in weak-coupling 
(^ ^ 1) region. It shows that the spin bath will weaken 








FIG. 2: (Color online) Concurrence evolution versus A for 
A^ = 801, J = 2, 7 = 1, D = when there is weak-coupling 
between qubits and spin bath, i.e., g — 0.05. 



the entanglement generation when the coupling between 
system and environment is adiabatically turned on. Es- 
pecially, the quantum phase transition of XY model at 
A = 1 will greatly enhance the decoherence, which coin- 
cides with the rapid decay of LE in an Ising model [l^ • 
It should be noted that we consider dynamics evolution 
of the system based on a finite-sized environment. To 
examine the effect of the chain length N on the quan- 
tum entanglement evolution, we plot Fig. [3l It shows 
that the longer the spin chain, the more serious the at- 
tenuation at QCP. The concurrence displays oscillatory 
decay of time for A/" = 401 and 801. As the length of 
the chain increases, the maximum value of concurrence 
decreases, and the revival of the concurrence disappears. 
The inset of Fig. [3] shows that the maximum concur- 
rence decays exponentially with the square root of chain 
length N . In addition, in weak-coupling region, for large 
A, the concurrence restores the sine function versus time. 
It seems that the effect of system-bath coupling is in- 
significant. The reason is that the spin bath is polarized 
along z axis in strong magnetic field. In this case, each 
spin is not entangled with the rest spins and the qubits 
(37I , [38| . In other word, the ferromagnetic bath can be 
more or less thought of as classical. Therefore, it has 
negligible effect on qubits, and the interaction between 
qubits plays a dominant role in entanglement generation 
again, as shown in Fig. [2l 

Then we investigate the effect of DM interaction at 
QCP of the spin bath. In Fig. |4] (a)-(c), we display 
the concurrence against time t under different D with 
7 equals to 0.8, 0.6 and 0.4, respectively. We can see 
the maximum concurrence is enhanced by increasing D. 
The bigger the magnitude of the less time needed 
to reach the maximum concurrence. The growth of en- 
tanglement can be interpreted via the fact that the DM 
interaction arouses the strong planar quantum fiuctua- 
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FIG. 3: (Color online) Concurrence versus time t is plotted tor 
different chain length N with A = 1, g=0.05, 7=1, J=2, D=0. 
From above to below, the lines correspond to N= 401, 801, 
1201, 1601 and 2001. The logarithm of maximum concurrence 
versus a/]V is shown in the inset. 




FIG. 4: (Color online) Concurrence as a function of time t 
under different D with (a) 7 — 0.8, (b) 7 = 0.6, and (c) 
7 — 0.4. We set other parameters as A = 1.0, J = 2, = 0.05, 
and N = 2001. 



tions. We would like to point out that the entanglement 
generation is suppressed at the QCP of spin bath as in 
Fig. [21 so the maximum concurrence can not reach unity. 
However, as the anisotropy parameter 7 decreases, the 
quantum fluctuations arising from the in-planar interac- 
tions become stronger [39], and then the two qubits in 
the presence of quantum fluctuations are more quantum 
correlated. In a word, the suppression of the entangle- 
ment at the QCP can be compensated both by increasing 
the DM interaction and by decreasing the anisotropy of 
the spin chain. 

We now turn to the entanglement dynamics beyond 
the weak-couphng region. Fig. [5] shows the concurrence 
as a function of time and the coupling g. When the 
coupling g increases, due to the decoherence of spin en- 
vironment, the concurrence damps dramatically. Thus 
in the strong-coupling region, no entanglement is gener- 
ated. However, it should be noted that in the vicinity of 
g = X = 1, the suppression of the entanglement genera- 
tion at the quantum critical point is released. At some 
time entanglement emerges suddenly and reaches a max- 
imum value in a short time. It can be easily concluded 
from Fig. [6] that the resonance peaks periodically appear 
at t = m7r/2, where m is an integer. We study the be- 
havior of the resonance peaks with respect to parameters 
D and 7. As shown in Fig. [6] (a), resonance peaks disap- 
pear quickly as 7 decreases, and there is no revival when 
7 is less than certain threshold. A similar phenomenon 
builds up as the parameter D increases, as shown in Fig. 
[6](b). This reveals that these resonance peaks are highly 
sensitive to both the anisotropy parameter and the DM 
interaction strength. 

In order to make a scrutiny into the resonance peaks, 
we rule out the effect of anisotropy and DM terms, i.e.. 




FIG. 5: (Color online) Concurrence as a function of time t 
and coupling parameter g with N = 801, where we set other 
parameters as A = 1.0, J = 2.0, 7=1, and D = 0. 



I) = 0, 7 = 1. Now we have Ai = 2 A, A2 =A3 = A, 
and A4 = 0. Under these parameters. Fig. [3 depicts 
the decoherence factors and concurrence against time in 
the resonance case ^ = A = 1. We can see that the 
decoherence factors |Fi2| and |Fi4| collapse to zero fleetly, 
while IF24I displays periodic revivals as time goes on. 
According to Eq. (|23|) . the concurrence is dominated by 
IF24I when IF12I and |Fi4| decay to zero, and exhibits the 
similar behavior as IF24I. 

To illustrate why there are robust peaks of concurrence 
for 7 = 1 and I) = 0, let us give a heuristical explanation 
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i^24(m7r/2)=n/e>o(-l)'^e^'^^^'^/^. The norm of 
F24(m7r/2) is unity. Based on such analysis, at t = m7r/2, 
we can approximatively set F12 = 0, F14 = 0, F24 = 1 in 
Eq. ([22]) , and thus we obtain C(m7r/2) = 0.5. For a tiny 
time deviation 6 from resonance time m7r/2, it arrives 
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FIG. 6: (Color online) Concurrence as a function of time t 
and (a) the anisotropy parameter 7, (b) the intensity of the 
DM interaction D. We set other parameters as N — 801, 
^ = A = 1.0, J 2.0 with (a) i:) = and (b) 7 = 1. 
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with A = X]/i.>o(-'-~^^^(^®/c^))- other words, the deco- 
herence factor F24(t) will exponentially decay in a short 
time after deviating from the resonance time. Therefore, 
at t 7^ m7r/2, replacing F12, F14 and F24 with zeros in 
Eq. ([22|) . we get = 0. It means that the entangle- 
ment generation is due to the presence of the momentum- 
independent spectrum structure. The flat band will also 
be destroyed when 7 is deviated from 1 or D is turned 
on, so the resonance peaks will be suppressed by varying 
7 or as shown in Fig. [6l These interesting features 
are also confirmed in our numerical study. With regard 
to the parameters of 7 = 1, = and ^ = A, we obtain 



A: 



/c>0 



2A2 sin^ k 



2Acos /c + 1 



2 ' 



A < 1, 
A > 1. 



(27) 
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(Color online) Docoherence factor |F/xi/| and concur- 
We set other parameters as 
i:) = and 7 = 1. 



FIG. 7 

rence as a function of time t. 
AT = 801, ^ = A = 1.0, J = 2.0, 



as follows. By Eq. (|2T]) . we have 



k>0 



(24 



Such form of expression has been studied carefully in 
Refs. [13, EqI, 113 • 1^12(^)1 shows a Gaussian decay with 
time for A = 1. Interestingly, the decoherence factor 
Fii{t) = 

x{G\e^^E t^-iH^ ^\G)x is similar to the decoherence fac- 
tor in the spin-echo experiment in Ref. [41], and also the 
decay of |Fi4(t)| is Gaussian [27], as shown in Fig. [Tja). 

A remarkable difference arises from F2A{t). From Eq. 
(|2T]) . we have 



Thus, the width of spikes of IF24I in Fig. Wljo) is inversely 
proportional to V^A when A < 1, while it is inversely 
proportional to VN when A > 1. Since the concurrence 
has the same dependence as F24, the resonance peaks at 
QCP will get sharper with the increase of chain length 
N. The peaks are robust even at the thermodynamic 
limit. 

We want to emphasize that the behavior of 1^24(^)1 is 
universal for any ^ = A, while the Gaussian decay for 
^12(14) (^) is only shown around QCP. When g = X ^ 
1, from the angle for Bogoliubov transformation in Eq. 
(p3|). we have 

^0, ^^«o, e, 



0, 



0, 



k, 

k. 



-TT < fc < 0, 
< fc < TT. 



(28) 
(29) 
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Consequently, 

Fi2{t) « n 

k>0 

Fuit) « II' 

k>0 



COS' 



sm' 



(30) 



(31) 



Due to the condition that g = \ 
nian has an important feature: the spectra Vt 
the modes are independent of momentum 
= 2. Consequently, at t = m7r/2 (m =0, 



all 



such Hamilto- 
spectra 1^^^ of 
i.e., 
1, 2, 



(25) Apparently, \Fi2(t)\ ^ 1, Fi^{t) / F24{t) - Fi2{t). It is 
notable that the form of Eq. (|3T]) is similar with Eq. 
([25|) . Compared with \F24{t)\ analyzed above, |Fi2(^)| 
should behave analogically, and spikes also spring up 
at t = m7r/2. The situation is contrary to the case of 
^ = A = 1. At t m7r/2, all the decoherence factors 
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IV. CONCLUSIONS 
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FIG. 8: (Color online) Concurrence as a function of time t for 
(a) N=201, (b) N=801, (c) N=1601 with ^=A=100, J = 2, 
7 = 1, and i:> = 0; (d) ^=A=10, (e) g=A=50, (f) g=A=500 
with N = 201, J = 2, 7 = 1, and i:> = 0. 



are zeros except Fi2(^), so the concurrence is dominated 
by |i^i2(^)|, and thus the concurrence forms a flat plat- 
form structure, which remains at = 0.5. As for 
t = m7r/2, the concurrence shows exotic patterns, which 
are sensitive to the chain length N and the amplitude of 
coupling g. Peaks and valleys emit from the platform at 
t = m7r/2, where the system can reach maximum or zero 
concurrence state periodically, as depicted in FiglHl Ac- 
cording to the analysis above, we know that the width of 
peaks and valleys is inversely proportional to It is 
observed that the width of peaks and valleys decreases as 
the chain length N increase (see Fig. [8](a)-(c)), and it is 
irrelevant to the couphng strength g (see Fig. [8](d)-(f)). 
For the region out of scope of QCP and strong coupling 
with ^ = A, there exists competitive effect among these 
three decoherence factors. 



The investigation of decoherence of spin environment is 
an attractive topic in the entanglement dynamical pro- 
cess. In the paper, we have studied the entanglement 
dynamics in term of concurrence of two interacting spin 
qubits coupled to a general XY spin chain with DM 
interaction. When there is no coupling between two- 
qubit system and environment, the interaction between 
the qubits will generate entanglement as a sine function 
versus time. In the weak-coupling region, the entan- 
glement generation is suppressed remarkably near the 
QCP of the spin bath. However, this suppression can 
be compensated both by increasing the DM interaction 
or by decreasing the anisotropy of the spin chain. Be- 
yond the weak-coupling region, the system-bath coupling 
will make the entanglement in the system decay rapidly 
except for ^ = A, where periodic resonance peaks of con- 
currence appear. By analyzing the decoherence factors, 
we attribute the appearance of resonance peaks to the 
emergence of the constant energy flat band. In such 
circumstances, increasing DM interaction or decreasing 
anisotropy will spoil the flatness of the band, and then 
the resonance peaks collapse. The effects of 7 and D 
here seem contrary to case of weak coupling. The seem- 
ingly contradictory phenomena are originated from dif- 
ferent mechanisms. The investigation of such model, in 
our opinion, may lead to several interesting phenomena 
in understanding decoherence in quantum open system. 
We hope that the simple analytical model described in 
this paper will assist us to gain further insight in QIP. 
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